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Abstract 

In i-theory a typical layer of a hierarchical architecture consists of HW modules pooling 
the dot products of the inputs to the layer with the transformations of a few templates under 
a group. Such layers include as special cases the convolutional layers of Deep Convolutional 
Networks (DCNs) as well as the non-convolutional layers (when the group contains only 
the identity). Rectifying nonlinearities - which are used by present-day DCNs - are one 
of the several nonlinearities admitted by i-theory for the HW module. We discuss here 
the equivalence between group averages of linear combinations of rectifying nonlinearities 
and an associated kernel. This property implies that present-day DCNs can be exactly 
equivalent to a hierarchy of kernel machines with pooling and non-pooling layers. Finally, 
we describe a conjecture for theoretically understanding hierarchies of such modules. A 
main consequence of the conjecture is that hierarchies of trained HW modules minimize 
memory requirements while computing a selective and invariant representation. 


1 Introduction 

The architectures now called Deep Learning Convolutional networks appeared with the name 
of convolutional neural networks in the 1990s - though the supervised optimization techniques 
used for training them have changed somewhat in the meantime. Such architectures have a 
history that goes back to the original Hubei and Wiesel proposal of a hierarchical architecture 
for the visual ventral cortex iterating in different layers the motif of simple and complex cells 
in VI. This idea led to a series of quantitative, convolutional cortical models from Fukushima 
(HI) to HMAX (Riesenhuber and Poggio, |2|). In later versions (Serre et al., [3j) such models 
of primate visual cortex have achieved object recognition performance at the level of rapid 
human categorization. More recently, deep learning convolutional networks trained with very 
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large labeled datasets (Russakovsky et al. |4|, Google |5j, Zeiler and Fergus 0) have achieved 
impressive performance in vision and speech classification tasks. The performance of these 
systems is ironically matched by our present ignorance of why they work as well as they do. 
Models are not enough. A theory is required for a satisfactory explanation and for showing the 
way towards further progress. This brief note outlines a framework towards the goal of a full 
theory. 

Its organization is as follows. We first discuss how i-theory applies to existing DLCNs. We 
then show that linear combinations of rectification stages can be equivalent to kernels. Deep 
Learning Convolutional Networks can be similar to hierarchies of HBFs (|7|). 

2 DCNs are hierarchies of kernel machines 

In this section, we review the basic computational units composing deep learning architectures 
of the convolution type. Then we establish some of their mathematical properties by using i- 
theory as described in f8j(9, 10]]. 

2.1 DCNs and i-theory 

The class of learning algorithms called deep learning, and in particular convolutional net¬ 
works, are based on a basic operation in multiple layers. We describe it using the notation of 
i-theory. 

The operation is the inner product of an input with another point called a template (or a 
filter, or a kernel), followed by a non linearity, followed by a group average. The output of the 
first two steps can be seen to roughly correspond to the neural response of a so called simple 
cell |ll 1| |. The collection of inner products of a given input with a template and its transfor¬ 
mations in i-theory corresponds to a so called convolutional layer in DCNs. More precisely, 
given a template t and its transformations gt, here g 6 Q is a finite set of transformations (in 
DCNs the only transformations presently used are translations), we have that each input x is 
mapped to ( x,gt ), g^Q- The values are hence processed via a non linear activation function, 
e.g. a sigmoid (1 + e -5 ) -1 , or a rectifier |s + b | + = max{-F,s} for s, b e 1R (the rectifier nonlinearity 
was called ramp by Breiman |ll21| ). In summary, the first operation unit can be described, for 
example by 

x\-+\(x,gt) + b\+. 

The last step, often called pooling, aggregates in a single output the values of the different 
inner products previously computed that correspond to transformations of the same template, 
for example via a sum 

^ 1 < x,gt ) + b\ + , t eT,b elR (1) 

g 

or a max operation 

max|(x, gt) + b\ + , teT,belR (2) 

g 

This corresponds to the neural response of a so called complex cell in [111] , 

2.2 A HW module is a kernel machine 

It is trivial that almost any (positive or negative defined) nonlinearity after the dot product in a 
network yields a kernel. Consider a 3-layers network with the first layer being the input layer 
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x. Unit i in the second layer (comprising N units) computes |(f;,x) + b t :| + = < pi(x),i = !,■■■ ,N. 
Thus each unit in the third layer performing a dot product of the vector of activities from 
the second layer with weights <pj(y) (from layer two to three) computes K(x,y ) = Jfj <Pj( x )<Pj(y) 
which is guaranteed to be well defined because the sum is finite. Key steps of the formal proof, 
which holds also in the general case, can be found in Rosasco and Poggio, |[101 | and references 
there; see also Appendix [5TT] and Appendix|5.2|for an example. 


A different argument shows that it is "easy" to obtain kernels in layered networks. Assume 
that inputs x e 1R ,! as well as the weights t are normalized, that is x, t e S n where S n is the unit 
sphere. 

In this case dot products are radial functions (since (x, t ) = ^(2-(|x-t| 2 ))) of r 2 . The kernel r 2 can 
be shaped by linear combinations (with bias) of rectifier units (linear combinations of ramps 
can generate sigmoidal units, linear combinations of which are equivalent to quasi-Gaussian, 
"triangular" functions )[] 

Thus for normalized inputs dot products with nonlinearities can easily be equivalent to 
radial kernels. Pooling before a similarity operation (thus pooling at layer n in a DCN before 
dot products in layer n + 1) maintains the kernel structure (since K(x,x') - J dgj dg'K(gx, g'x') 
is a kernel if K is a kernel). Thus DCNs with normalized inputs are hierarchies of radial kernel 
machines, also called Radial Basis Functions (RBFs). 

Kernels induced by linear combinations of features such as 


Y^\(x,gt) + b\ + , t sT,b 6 1R 


(3) 


are selective and invariant (if G is compact, see Appendix 5.1 >. The max operation max^ | (x,gt)+ 
b\ + has a different form and does not satisfy the condition of the theorems in the Appendix. On 
the other hand the pooling defined in terms of the following soft-max operation (which ap¬ 
proximates the max for "large" n) 


L 


« x,gt)) n 


LAi + ix.g't))”- 1 ’ 


teT 


(4) 


induces a kernel that satisfies the conditions of the theorems summarized in Appendix 5.1 


It is well known that such an operation can be implemented by simple circuits of the lateral 
inhibition type (see HU). On the other hand the pooling in Appendix |5.3| (see [IT5l| l does not 
correspond in general to a kernel, does not satisfy the conditions of Appendix |5. 1 1 and is not 
guaranteed to be selective. 

Since weights and "centers" of the RBFs are learned in a supervised way, the kernel ma¬ 
chines should be more properly called HyperBF, see Appendix |5.4| 


2.3 Summary: every layer of a DCN is a kernel machine 

Layers of a Deep Convolutional Network using linear rectifiers (ramps) can be described as 

I ( x >gt) + b\ + , t eT,b elR (5) 

g 

Notice that in one dimension the kernel \x - y | can be written in terms of ramp functions as \x - y\ = \x - y| + +1 - 
(x-y)|+. See Figure[l]and (TTj. 
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Figure 1: A “Gaussian” of one variable can be -written as the linear combination of ramps (e.g. 
rectifiers): a) a sigmoid-like function can be written as linear combinations of ramps b) linear com¬ 
binations of sigmoids give gaussian-like triangular functions. 


where the range of pooling (y ? ) may be degenerate (no pooling) in which case the layer is not 
a convolutional layer. 

Such a layer corresponds to the kernel 

K(x,x') = J dg dg' K 0 (x,g,g'). (6) 

with 


K o( x ’g’g') = J dtdb | (gt,x) + b\ + \(g’t,x’) + b\ + . 

An alternative pooling to [5]is provided by softmax pooling 


(7) 


L^j^ g ,(i + ( x ,g't)) n - 1 


( 8 ) 


Present-day DCNs seem to use equation[8]for the pooling layers and equation [5]for the non- 
convolutional layers. The latter is the degenerate case of equation [8] (when G contains only the 
identity element). 


3 A conjecture on what a hierarchy does 

The previous sections describe an extension of i-theory that can be applied exactly to any layer 
of a DLCN and any of the nonlinearities that have been used: pooling, linear rectifier, sigmoidal 
units. In this section we suggest a framework for understanding hierarchies of such modules, 
which we hope may lead to new formal or empirical results. 
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Figure 2: A hierarchical, supervised architecture built from eHW-modules. Each red circle represents 
the signature vector computed by the associated module (the outputs of complex cells) and double 
arrows represent its receptive fields - the part of the (neural) image visible to the module (for trans¬ 
lations this is also the pooling range). The “image” is at level 0, at the bottom. 

3.1 DLCNs and HVQ 

We consider here an HW module with pooling (no pooling corresponds to 1 pixel stride con¬ 
volutional layers in DLN). Under the assumption of normalized inputs, this is equivalent to 
a HBF module with Gaussian-like radial kernel and "movable" centers (we assume standard 
Euclidean distance, see Poggio, Girosi.p~6|; Jones, Girosi and Poggio,p~3|). 

Notice that one-hidden-layer HBF can be much more efficient in storage (e.g. bits used for 
all the centers) than classical RBF because of the smaller number of centers (HBFs are similar 
to a multidimensional free-knots spline whereas RBFs correspond to classical spline). 

The next step in the argument is the observation that a network of radial Gaussian-like 
units become in the limit of a —> 0 a look-up table with entries corresponding to the centers. 
The network can be described in terms of soft Vector Quantization (VQ) (see section 6.3 in Pog¬ 
gio and Girosi, {7|). Notice that hierarchical VQ (dubbed HVQ) can be even more efficient than 
VQ in terms of storage requirements (see e.g. [1171)1. This suggests that a hierarchy of HBF layers 
may be similar (depending on which weights are determined by learning) to HVQ. Note that 
compression is achieved when parts can be reused in higher level layers as in convolutional networks. 
Notice that the center of one unit at level n of the "convolutional" hierarchy of Figure [2] is a 
combinations of parts provided by each of the lower units feeding in it. This may even happen 
without convolution and pooling as shown in the following extreme example. 

Example Consider the case of kernels that are in the limit delta-like functions (such as Gaus¬ 
sian with very small variance). Suppose as in Figure[3]that there are four possible quantizations 
of the input x: x 4 ,x 2 ,x 3 ,x 4 . One hidden layer would consist of four units <5(x-x ; ),z = l,--- ,4. 
But suppose that the vectors x lr x 2 ,x 3 ,x 4 can be decomposed in terms of two smaller parts or 
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(A, A) ( A,B ) (B,A) (S,B) 


X 1 X 2 X 3 X4 




Figure 3: See text, Example in section 
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features x’ and x", e.g. X\ = x' ®x”, x 2 = x’ ®x’, x 3 = x" ®x” and X 4 = x” ®x'. Then a two layer 
network could have two types of units in the first layer S(x - x') and 5(x - x”); in the second 
layer four units will detect the conjunctions of x' and x" corresponding to Xi,x 2 ,x 3 ,X 4 . The 
memory requirements will go from 4 N to 2N/2 + 8 where N is the length of the quantized 
vectors; the latter is much smaller for large N. Memory compression for HVQ vs VQ - that is 
for multilayer networks vs one-layer networks - increases with the number of (reusable) parts. 
Thus for problems that are compositional, such as text and images, hierarchical architectures of 
HBF modules minimize memory requirements. 

Classical theorems (see refrences in [18.Il9lj show that one hidden layer networks can ap¬ 
proximate arbitrarily well rather general classes of functions. A possible advantage of multi¬ 
layer vs one-layer networks that emerges from the analysis of this paper is memory efficiency 
which can be critical for large data sets and is related to generalization rates. 

4 Remarks 

• Throughout this note, we discuss the potential properties of multilayer networks, that is 
the properties they have with the "appropriate" sets of weights when supervised training 
is involved. The assumption is therefore that greedy SGD using very large sets of labeled 
data, can find the "appropriate sets of sets of weights". 

• Recently, several authors have expressed surprise when observing that the last hidden 
unit layer contains information about tasks different from the training one (e.g. [l20l|l. 
This is in fact to be expected. The last layer of HBF is rather independent of the training 
target and mostly depends on the input part of the training set (see theory and gradient 
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descent equations in Poggio and Girosi, [21 for the one-hidden layer case). This is ex¬ 
actly true for one-hidden layer RBF networks and holds approximatively for HBFs. The 
weights from the last hidden layer to the output are instead task/target dependent. 

• The result that linear combinations of rectifiers can be equivalent to kernels is robust in 
the sense that it is true for several different nonlinearities such as rectifiers, sigmoids etc. 
Ramps (e.g. rectifiers) are the most basic ones. Such robustness is especially attractive for 
neuroscience. 
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5 Appendices 

5.1 HW modules are equivalent to kernel machines (a summary of the results in 

m) 

In the following we summarize the key passages of [Il0| in proving that HW modules are kernel 
machines: 

1. The feature map 

c[)(x, t,b ) = \(t,x) + b\ + 

(that can be associated to the output of a simple cell, or the basic computational unit of 
a deep learning architecture) can also be seen as a kernel in itself. The kernel can be a 
universal kernel. In fact, under the hypothesis of normalization of the vectors x, t we have 
that 2(| 1 ~(t,x)\ + + \(t,x)~ 1| + ) = 2|1 -(t,x)\ = \\x- t||| which is a universal kernel (see also 
th 17 of [Ell). 

The feature <fi leads to a kernel 

K 0 (x,x') = cp T (x)cp(x') - db dt \(t,x) + b\ + \(t,x') + b\ + 

which is a universal kernel being a kernel mean embedding (w.r.t. t, b, see |I22I| ) of the a 
product of universal kernels. 

2. If we explicitly introduce a group of transformations acting on the feature map input i.e. 
< p(x,g , t, b) = \(gt,x) + b | + the associated kernel can be written as 

dgdg' J dt db\(gt,x) + b\ + \(g't,x') + b\ + = J dg dg' K 0 (x,g,g'). 

K(x,x') is the group average of K 0 (see [I23lj ) and can be seen as the mean kernel embed¬ 
ding of K 0 (w.r.t. g,g', see f22]). 

3. The kernel K is invariant and, if G is compact, selective i.e. 

K(x,x') = 1 o x ~ x . 

The invariance follows from the fact that any G-group average function is invariant to G 
transformations. Selectivity follows from the fact that K a universal kernel being a kernel 
mean embedding of K 0 which is a universal kernel (see ||22|| ). 

Remark 1. If the distribution of the templates t follows a gaussian lau> the kernel K 0 , with an oppor¬ 
tune change of variable, can be seen as a particular case of the nth order arc-cosine kernel in H24\l for 
n = 1. 



5.2 An example of an explicit calculation for an inner product kernel 

Here we note how a simple similarity measure between functions of the form in eq. § cor¬ 
respond to a kernel in the case when the inner product between two HW module outputs at 
the first layer, say p(I),p(I'), is calculated using a step function nonlinearity. Note first that 
a heaviside step function can be approximated by a sigmoid like function derived by a linear 
combination of rectifiers of the form: 

H(x)~a( |x| + -|x--| + ) 
a 
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for very large values of a. With this specific choice of the nonlinearity we have that the inner 
product of the HW modules outputs (for fixed t) at the first layer is given by: 

<p(I),/*(!')> = \\ d S d s\\ d bH(b - (I,gt))H(b - < I'.g't »). 

with 

Fbi 1 ) = J d g H (b-(I,gt)). 

Assuming that the scalar products, (I,gt),(I',g't), range in the interval [-p,p] a direct compu¬ 
tation of the integral above by parts shows that, x = ( I,gt),x’ = (I',gt): 


r 


dbH(b - x)H(b - x') 


= H(b - x)((b - x')H(b - x') - (x-x')H(x-x' 


: p - —(x + x’ + \x-x’\) 



and being 


max{x,x'} = max|x- —(x + x),x - —(x + x')j + —(x + x') 

= max {\( x '- x )’\( x ~ x ')} + \( x + x ') 

= max (--(x- x'), X -{x- x')} + \{x + x') 


= + 




+ -(x + x) 
2 ’ 


- —(x + X + \x — x'\). 


we showed that 

K((I\t), (I, t)) = C- max«r, t), (I, t)) 

which defines a kernel. If we include the pooling over the transformations and templates we 
have 

£(/,/') = P-J dA{t)dgdg' max((I',gt),(I,g't)). (9) 

where A(f) is a probability measure on the templates t. K is again a kernel. 

A similar calculation can be repeated at the successive layer leading to kernels of kernels struc¬ 
ture. 


5.3 Mex 

Mex is a generalization of the pooling function. From [II5j eq. 1 it is defined as: 


Mex, 


(k,U) 


^ lo g(^L ex p ( ^)) 

2 = 1 


( 10 ) 
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We have 


Mgx ({c.-W) T- > Maxj(ci) 

C, —>00 

Mex ((c;U) Meani(cj) 

Mex ({cM) ~ - > Mirii(ci). 

c ,—>—oo 

We can also choose values of <5; in between the ones above, the interpretation is less obvious. 
The Mex pooling does not define a kernel since is not positive definite in general (see also Th 1 
in H3) 


5.4 Hyper Basis Functions: minimizing memory in Radial Basis Function net¬ 
works 


We summarize here an old extension by Poggio and Girosi |l25l| of the classical kernel networks 
called Radial Basis Functions (RBF). In summary (but see the paper) they extended the theory 
by defining a general form of these networks which they call Hyper Basis Functions. They 
have two sets of modifiable parameters: moving centers and adjustable norm-weights. Moving the 
centers is equivalent to task-dependent clustering and changing the norm weights is equivalent 
to task-dependent dimensionality reduction. 

A classical RBF has the form 


N 


fix) = ^C;K(||x-Xi|| 2 ), 
1=1 


( 11 ) 


which is a sum of radial functions, each with its center x ( on a distinct data point. Thus the 
number of radial functions, and corresponding centers, is the same as the number of examples. 
Eq. (|11[> is a minimizer solution of 


N 

H[/] = ^-/(x i )) 2 + A||F/|| 2 A 6 IR + 

i =1 


( 12 ) 


where P is a constrain operator (usually a differential operator). 

HBF extend RBF in two directions: 

1. The computation of a solution of the form <du has a complexity (number of radial func¬ 
tions) that is independent of the dimensionality of the input space but is on the order 
of the dimensionality of the training set (number of examples), which can be very high. 
Poggio and Girosi showed how to justify an approximation of equation ( |11| in which the 
number of centers is much smaller than the number of examples and the positions of the 
centers are modified during learning. The key idea is to consider a specific form of an 
approximation to the solution of the standard regularization problem. 

2. Moving centers are equivalent to the free knots of nonlinear splines. In the context of 
networks they were first suggested as a potentially useful heuristics by Broomhead and 
Lowe |26j and used by Moody and Darken [27]. 

Poggio and Girosi called Hyper Basis Functions, in short HyperBFs, the most general form of 
regularization networks based on these extensions plus the use of a weighted norm. 
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5.4.1 Moving Centers 

The solution given by standard regularization theory to the approximation problem can be 
very expensive in computational terms when the number of examples is very high. The com¬ 
putation of the coefficients of the expansion can become then a very time consuming operation: 
its complexity grows polynomially with N, (roughly as N 3 ) since an N x N matrix has to be in¬ 
verted. In addition, the probability of ill-conditioning is higher for larger and larger matrices 
(it grows like N 3 for a N x N uniformly distributed random matrix) [l28l |. The way suggested 
by Poggio and Girosi to reduce the complexity of the problem is as follows. While the exact 
regularization solution is equivalent to generalized splines with fixed knots, the approximated 
solution is equivalent to generalized splines with/ree knots. 

A standard technique, sometimes known as Galerkin's method, that has been used to find 
approximate solutions of variational problems, is to expand the solution on a finite basis. The 
approximated solution /*(x) has then the following form: 

n 

r(x) = j^ c ict> i (x) (i3) 

i=l 

where {<piYl =1 is a set of linearly independent functions [l29lj . The coefficients c,- are usually 
found according to some rule that guarantees a minimum deviation from the true solution. A 
natural approximation to the exact solution will be then of the form: 

n 

/*(x) = ^c a G(x;t a ) (14) 

a-1 

where the parameters t a , that we call "centers", and the coefficients c a are unknown, and are 
in general fewer than the data points (n < N). This form of solution has the desirable property 
of being an universal approximator for continuous functions (7| and to be the only choice that 
guarantees that in the case of n - N and {t a }” =1 = {x ; }” =| the correct solution (of equation 
) is consistently recovered. We will see later how to find the unknown parameters of this 
expansion. 


5.4.2 How to learn centers' positions 

Suppose that we look for an approximated solution of the regularization problem of the form 

n 

f*(x)=J^ c aG(\\x-ta\\ 2 ) (15) 

a=1 


We now have the problem of finding the n coefficients c a , the dx n coordinates of the centers 
t a . We can use the natural definition of optimality given by the functional H. We then impose 
the condition that the set {c a ,t a \a = l,...,n} must be such that they minimizes H[f*], and the 
following equations must be satisfied: 


dH[H 

dc a 


= 0 , 


dHlT 

dt„ 


0, a = 1, 


n . 


Gradient-descent is probably the simplest approach for attempting to find the solution to this 
problem, though, of course, it is not guaranteed to converge. Several other iterative methods, 
such as versions of conjugate gradient and simulated annealing |l30f may be more efficient 
than gradient descent and should be used in practice. Since the function H[f*] to minimize is 
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in general non-convex, a stochastic term in the gradient descent equations may be advisable 
to avoid local minima. In the stochastic gradient descent method the values of c a , t a and M 
that minimize H[f*] are regarded as the coordinates of the stable fixed point of the following 
stochastic dynamical system: 


- 0 )- 


dH[f* 

dc„ 


+ rja(t)> a = l,...,n 


dH[f* 1 

= — + F a ( t )> a = l,...,n 

where t] a (t), f* a (t) are white noise of zero mean and to is a parameter determining the micro¬ 
scopic timescale of the problem and is related to the rate of convergence to the fixed point. 
Defining 

n 

\ = yi-f*(^) = Vi-J^ c aG(\\Xi-t a \\ 2 ) 
a -1 


we obtain 


N 

h[T] = h cX = Y} a,) 2 . 

z=1 


The important quantities - that can be used in more efficient schemes than gradient descent - 
are 


• for the c a 


• for the centers t a 


mn 

dc a 


N 

-2^A,G(||x i -t„|| 2 ) ; 

i =1 


mn 

dt a 


N 

= 4 c a ^A ; G'(||x ! -t a || 2 )(x ! - t a ) 

z=l 


(16) 

(17) 


Remarks 

1. Equation ( |16| has a simple interpretation: the correction is equal to the sum over the 
examples of the products between the error on that example and the "activity" of the 
"unit" that represents with its center that example. Notice that H[f*] is quadratic in the 
coefficients c a , and if the centers are kept fixed, it can be shown |7J that the optimal 
coefficients are given by 

c = (G t G + Ag)~ 1 G T y (18) 

where we have defined (y),- = (c) a = c a , (G) ;a = G(x,;t a ) and (g) a p = G(t a ;tg). If A is let 
go to zero, the matrix on the right side of equation ( p~8| ) converges to the pseudo-inverse of 
G |[31| | and if the Green's function is radial the approximation method of ||26l| is recovered. 

2. Equation ( |T7j l is similar to task-dependent clustering |7J. This can be best seen by as¬ 
suming that A, are constant: then the gradient descent updating rule makes the centers 
move as a function of the majority of the data, that is of the position of the clusters. In 
this case a technique similar to the k-means algorithm is recovered, [32] |27j. Equating 
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dH }/ - to zero we notice that the optimal centers i a satisfy the following set of nonlinear 

O'-a 

equations: 


t« = 


LiP?* 

L t P? 


a = l,...,n 


where P“ = A ( G'(||x ; - tj| 2 ). The optimal centers are then a weighted sum of the data 
points. The weight P" of the data point z for a given center t a is high if the interpolation 
error A, is high there and the radial basis function centered on that knot changes quickly 
in a neighborhood of the data point. This observation suggests faster update schemes, in 
which a suboptimal position of the centers is first found and then the c a are determined, 
similarly to the algorithm developed and tested successfully by Moody and Darken [l27l |. 


5.4.3 An algorithm 

It seems natural to try to find a reasonable initial value for the parameters c,t ff , to start the 
stochastic minimization process. In the absence of more specific prior information the follow¬ 
ing heuristics seems reasonable. 

• Set the number of centers and set the centers' positions to positions suggested by cluster 
analysis of the data (or more simply to a subset of the examples' positions). 

• Use matrix pseudo-inversion to find the c a . 

• Use the t a , and c a found so far as initial values for the stochastic gradient descent equa¬ 
tions. 

Experiments with movable centers and movable weights have been performed in the context of 
object recognition (Poggio and Edelman, [1331] ; Edelman and Poggio, f34| ) and approximation 
of multivariate functions. 


5.4.4 Remarks 

1. Equation ( ]17[ > is similar to a clustering process. 

2. In the case of N examples, n — N fixed centers, there are enough data to constrain the N 
coefficients c a to be found. Moving centers add another nd parameters (d is the number 
of input components). Thus the number of examples N must be sufficiently large to 
constrain adequately the free parameters - n d-dimensional centers, n coefficients c a . 
Thus 

N » n + nd. 


Acknowledgment 

This work was supported by the Center for Brains, Minds and Machines (CBMM), funded by 
NSF STC award CCF - 1231216. This work was also supported by A*STAR JCO VIP grant 
#1335h00098. Part of the work was done in Singapore at the Institute for Infocomm under 
REVIVE funding. TP thanks A*Star for its hospitality. 


13 




References 


[1] K. Fukushima, "Neocognitron: A self-organizing neural network model for a mecha¬ 
nism of pattern recognition unaffected by shift in position," Biological Cybernetics, vol. 36, 
pp. 193-202, Apr. 1980. 

[2] M. Riesenhuber and T. Poggio, "Hierarchical models of object recognition," Nature Neuro¬ 
science, vol. 3,11, 2000. 

[3] T. Serre, A. Oliva, and T. Poggio, "A feedforward architecture accounts for rapid catego¬ 
rization," Proceedings of the National Academy of Sciences of the United States of America, 
vol. 104, no. 15, pp. 6424-6429, 2007. 

[4] O. Russakovsky, J. Deng, H. Su, J. Krause, S. Satheesh, S. Ma, Z. Huang, A. Karpathy, 
A. Khosla, B. M., A. Berg, and F.-F. L., "Large scale visual recognition challenge," Ima- 
geNet, arXiv.l 409.0575, 2014. 

[5] C. Szegedy, W. Liu, Y. Jia, P. Sermanet, S. Reed, D. Anguelov, D. Erhan, V. Vanhoucke, and 
A. Rabinovich, "Going Deeper with Convolutions," ArXiv e-prints, Sept. 2014. 

[6] M. D. Zeiler and R. Fergus, "Visualizing and understanding convolutional networks.," 
CoRR, vol. abs/1311.2901, 2013. 

[7] T. Poggio and F. Girosi, "A theory of networks for approximation and learning," Labora¬ 
tory, Massachusetts Institute of Technology, vol. A.I. memo nll40, 1989. 

[8] F. Anselmi, J. Z. Leibo, L. Rosasco, J. Mutch, A. Tacchetti, and T. Poggio, "Unsuper¬ 
vised learning of invariant representations in hierarchical architectures," arXiv preprint 
arXiv.l 311.4158, 2013. 

[9] F. Anselmi, L. Rosasco, and T. Poggio, "On invariance and selectivity in representation 
learning," arXiv.l 503.05938 and CBMM memo n 29, 2015. 

[10] L. Rosasco and T. Poggio, "Convolutional layers build invariant and selective reproducing 
kernels," CBMM Memo, in preparation, 2015. 

[11] D. Hubei and T. Wiesel, "Receptive fields, binocular interaction and functional architec¬ 
ture in the cat's visual cortex," The Journal of Physiology, vol. 160, no. 1, p. 106, 1962. 

[12] L. Breiman, "Hinging hyperplanes for regression, classification, and function approxima¬ 
tion," Tech. Rep. 324, Department of Statistics University of California Berkeley, Califor¬ 
nia 94720, 1991. 

[13] F. Girosi, M. Jones, and T. Poggio, "Regularization theory and neural networks architec¬ 
tures," Neural Computation, vol. 7, pp. 219-269, 1995. 

[14] M. Kouh and T. Poggio, "A canonical neural circuit for cortical nonlinear operations," 
Neural computation, vol. 20, no. 6, pp. 1427-1451, 2008. 

[15] N. Cohen and A. Shashua, "Simnets: A generalization of convolutional networks," CoRR, 
vol. abs/1410.0781, 2014. 

[16] T. Poggio and F. Girosi, "Regularization algorithms for learning that are equivalent to 
multilayer networks," Science, vol. 247, pp. 978-982, 1990. 


14 



[17] J. Mihalik, "Hierarchical vector quantization, of images in transform domain.," ELEK- 
TROTECHN. CA5, 43, NO. 3. 92,94., 1992. 

[18] F. Girosi and T. Poggio, "Representation properties of networks: Kolmogorov's theorem is 
irrelevant," Neural Computation, vol. 1, no. 4, pp. 465-469, 1989. 

[19] F. Girosi and T. Poggio, "Networks and the best approximation property," vol. 63, pp. 169- 
176, 1990. 

[20] D. Yamins, H. Hong, C. Cadieu, E. Solomon, D. Seibert, and J. DiCarlo, "Performance- 
optimized hierarchical models predict neural responses in higher visual cortex," Proceed¬ 
ings of the National Academy of Sciences, 2014. 

[21] C. A. Micchelli, Y. Xu, and H. Zhang, "Universal kernels.," Journal of Machine Learning 
Research, vol. 6, pp. 2651-2667, 2006. 

[22] B. K. Sriperumbudur, A. Gretton, K. Fukumizu, B. Scholkopf, and G. R. G. Lanckriet, 
"Hilbert space embeddings and metrics on probability measures," Journal of Machine 
Learning Research, vol. 11, pp. 1517-1561,2010. 

[23] B. Haasdonk and H. Burkhardt, "Invariant kernel functions for pattern analysis and ma¬ 
chine learning," Mach. Learn., vol. 68, pp. 35-61, July 2007. 

[24] Y. Cho and L. K. Saul, "Kernel methods for deep learning," in Advances in Neural In¬ 
formation Processing Systems 22 (Y. Bengio, D. Schuurmans, J. Lafferty, C. Williams, and 
A. Culotta, eds.), pp. 342-350, Curran Associates, Inc., 2009. 

[25] T. Poggio and F. Girosi, "Extensions of a theory of networks for approximation and learn¬ 
ing: dimensionality reduction and clustering," Laboratory, Massachusetts Institute of Tech¬ 
nology, vol. A.I. memo n 1167, 1994. 

[26] D. S. Broomhead and D. Lowe, "Multivariable Functional Interpolation and Adaptive Net¬ 
works," Complex Systems 2, pp. 321-355, 1988. 

[27] J. Moody and C. J. Darken, "Fast Learning in Networks of Locally-Tuned Processing 
Units," Neural Computation, vol. 1, no. 2, pp. 281-294, 1989. 

[28] J. Demmel, "The geometry of ill-conditioning," /. Complexity, vol. 3, pp. 201-229, 1987. 

[29] S. Mikhlin, The problem of the minimum of a quadratic functional. San Francisco, CA: 
Holden-Day, 1965. 

[30] S. Kirkpatrick, C. Gelatt, and M. Vecchi, "Optimization by simulated annealing," Science, 
vol. 220, pp. 219-227, 1983. 

[31] A. E. Albert, Regression and the Moore-Penrose pseudoinverse. Mathematics in science and 
engineering, New York, London: Academic Press, 1972. 

[32] J. MacQueen, "Some methods of classification and analysis of multivariate observations," 
in Proc. 5th Berkeley Symposium on Math., Stat., and Prob. (L. LeCam and J. Neyman, eds.), 
p. 281, Berkeley, CA: U. California Press, 1967. 

[33] T. Poggio and S. Edelman, "A network that learns to recognize 3D objects," Nature, 
vol. 343, pp. 263-266, 1990. 


15 



[34] S. Edelman and T. Poggio, "Bringing the grandmother back into the picture: a memory- 
based view of object recognition," a.i. memo 1181, mitai, 1990. 


16 



